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1 Introduction 



Neutrino oscillations, if neutrinos are massive, were first proposed by Pontecorvo [ij and by Maki, 
et al., m, respectively, in the late 1950's and early f960's. It is to date the only experimental 
indication of physics beyond the standard model (SM) that may provide a starting point to further 
expand the horizon of our most basic knowledge of nature. Historically neutrinos have played 
mostly passive roles in particle physics research. However, the observation of neutrino oscillations 
by the Super-Kamiokande experiment |J, which are corroborated by various other experiments, 
changed significantly the role of neutrinos and it has very profound implications for particle physics, 
astrophysics and cosmology. Neutrinos are the end product of the decays of almost all elementary 
particles, observed or proposed. They fill the universe as a relic background radiation similar to 
the photon as a finger print of the history of the universe. Now they are taking the center stage in 
particle physics research and their heightened fascination by physicists has just began. 

Since it was proposed more than thirty years ago, the standard model of electroweak and 
strong interactions (SM) has been firmly established as the fundamental theoretical framework of 
elementary particles. Its applications straddles eleven orders of magnitude in energy range, from 
the atomic parity violation of 1 eV to the W, Z and top quark physics of hundreds of GeV. All 
high energy phenomenologies that have been scrutinized by experiments are found to agree with 
the predictions of the SM. The only missing piece is the Higgs particle. Even for that, a tentative 
gleam of its existence from the last LEP data 0] indicated that it may be what is expected from the 
SM. In reaching the present status of the SM, the IHEP's BEPC program has made the important 
contribution in establishing the lepton universality which is a corner stone of the SM. In this 
"business as usual" state of affairs, evidence for neutrino oscillation and hence massive neutrinos 
gives light to the possibility of a new point of departure for post SM physics. 

From the available neutrino oscillation data, we can already see some rich features of the 
neutrino sector. In analogy to the mass spectra of the quarks and charged leptons, there is a 
hierarchical structure in the mass square differences of the neutrinos. But in a departure from the 
feature of the quark sector, the mixing angles of the neutrinos are large, some probably maximal. 
Another intriguing aspect of the neutrino oscillations is that they are manifested macroscopically 
at terrestrial and solar distances. The physics of neutrino oscillation can be treated effectively by 
simple quantum mechanics. When finally confirmed, it will be an unique and elegant illustration 
of the effect of quantum interference at macroscopic scales. 

Although the existing data offer strong indication of neutrino oscillations, the oscillation pa- 
rameters have not been determined with sufficient accuracy. The unique signature of the flavor 
transmutation, i.e., the appearance of a flavor different from the original one, has not been con- 
vincingly observed. In the ongoing and next generation neutrino oscillation experiments under 
construction, some of the parameters will be probed with greater accuracy. But some other pa- 
rameters may not yet be accessible. Hence more experiments designed to look for the missing 
information and to probe the known parameters with even greater accuracy are necessary. 

One possible implication of the macroscopic manifestation of the neutrino oscillation is that 
the effect of the lepton flavor mixing may be very small in microscopic distances. Hence there are 
possibly very little observable effects to the conventional high energy physics parameters measured 
at microscopic distances, such as the charge lepton flavor changing neutral current ||. This pro- 
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vides additional motivation for pursuing refined neutrino oscillation experiments so that the full 
implication of neutrino flavor mixing can be thoroughly investigated and understood, although the 
search for the effects at microscopic distances should be pursued. 

This report summarizes the preliminary results of our study on a very long baseline (LBL) 
neutrino oscillation experiment, which uses the intensive conventional neutrino beam from the 
High Intensity Proton Accelerator (HIPA) Q, Japan, or the beam of a possible neutrino factory, 
delivered to a detector located in Beijing, China, tentatively called the Beijing Astrophysics and 
Neutrino Detector (BAND). HIPA to BAND (H2B) has a distance of 2100 km, which is eight times 
longer than the currently online LBL experiment, K2K, and three times longer than all of the 
three LBL experiments in construction: MINOS, ICARUS and OPERA. This very long baseline 
program has the distance to neutrino energy optimal for the atmospheric neutrino oscillation scale. 
It enables us to investigate some of the parameters that are not easily accessible by the above 
mentioned on-going and approved experiments. For a neutrino factory located in North America 
or Europe the distance to BAND is even much longer. 

By the time this experiment becomes online some of the basic parameters related to the solar 
and atmospheric neutrino oscillations will be accurately determined. BAND will be in position to 
carry out a refined investigation of the neutrino oscillation parameters, such as the mixing param- 
eters sin 2 2#i3, the matter effect and the sign of the dominant neutrino mass-squared difference. 
Additional important avenues of investigation include the effect of CP violation, some of the as- 
trophysics measurements, the neutral current reactions that have scanty experimental information 
available to date, etc. The recent evidence of direct CP violation in the hadron sector |7|] (and 
therefore evidence of non-vanishing phase in the CKM matrix) suggests that a non-vanishing phase 
angle may also exist in the lepton sector. Furthermore, the large mixing angles of the neutrinos 
makes it plausible to expect that the CP phase angle is likely large in the lepton sector. 

This report is motivated by the discussions took placed at several workshops held at IHEP 
during the past two years. Several Physicists from Japan and USA joined their Chinese colleagues 
in the workshop held on May 25, 2000 ||. As an outcome, an agreement was made to form two 
study groups, both consisting of experimentalists and theorists, to make independent studies of the 
physics potential of H2B. One group, coordinated by Dr. Kaoru Hagiwara, is made of physicists 
from KEK and the University of Kyoto and the other group, which is responsible for the present 
study report, includes physicists from China and USA. We believe our study so far and that of the 
KEK-Kyoto group || show that the physics of H2B is compelling and merits serious considerations 
for this very long baseline neutrino oscillation experiment. Further studies using more realistic 
neutrino beam parameters with detailed detector design are called for. 

We envisage that the experimental programs outlined in this report will take a significant lead 
time to be fully implemented. The long lead time required is partly due to the R&D efforts of 
BAND but mainly constrained by the timeline of the neutrino beam, such as that from HIPA or a 
neutrino factory. We should also note that the long lead time is characteristic of modern large high 
energy physics experimental programs. Anticipating the long lead time, BAND will be constructed 
in the flexible modular approach and designed to do astrophysics experiment in the first stage of 
the program. In the intervening time, BAND will be thoroughly tested, enlarged and improved to 
be ready for the exciting physics of a LBL neutrino experiment. 

In Sec. 2 we summarize the current experimental status of neutrino oscillations and the expec- 
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tations of the new generation of experiments online and in construction. In Sec. 3 we summarize 
the theory of neutrino oscillations and list some of useful formulae. Section 4 outlines the funda- 
mentals of the long baseline neutrino oscillation experiment and presents the result of a preliminary 
study of the physics expectations of H2B. Sec. 5 presents a preliminary design of a possible BAND 
detector as a concrete description of the physics results reachable at H2B. 

A number of topics which are not covered by the present report will be the focus of the future 
studies, which involve a refined study of possible detector(s), study of other physics opportunities 
including those listed in Sec. 5.3, oscillation physics at even long oscillation distances that may 
involve BAND, and the comparison of H2B with other LBL experiments in construction. Topics 
such as the near detector, the neutrino beam design, and the choice of a far detector are expected 
to be the joint effort of a full collaboration with the participation of many physicists world wide, 
including the present study group, the critical participation from Japan, and from institutions of 
many other regions and countries. 
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2 Current experimental status of neutrino oscillations 



In the following we review briefly the present experiments and the relevant data. For more details 
we refer to review articles available in the literature; we list a few in Ref . |10| . We do not claim to be 
exhaustive and apologize for the omission of any articles and experimental data. For an exhaustive 
list of neutrino experiments see the neutrino web-site in Ref. [11]. Since the interpretation of 
the data are mostly in the two-flavor scenario, we will first summarize the formalism of two-flavor 
oscillations in vacuum. 



2.1 Two-flavor neutrino oscillations 

We begin in the lepton flavor framework in which the charged lepton mass matrix is diagonalized. 
Define the two-flavor eigenstates v a and u», and their mass eigenstates, v\ and v<i of masses m\ and 
m2, respectively. The flavor states and mass eigenstates are related by a mixing matrix U, which 
is an orthogonal transformation in two dimensions: 




where 9 is the mixing angle. The states are orthonormalized within their own spaces, i,e., 

< Vj\v k >= S jk , j,k = a,/3; or 1,2. (2) 
It should be noted that in an experiment, neutrinos are always produced as flavor eigenstates. 



The time evolution of a flavor state can be simply expressed in terms of the time evolution of 
the mass eigenstates which enter into the flavor state at t = 0, 

Suppose the neutrino flavor state u a is produced, then at time t we have 

u a (t) = (cos 2 9e' iElt + sin 2 9e~ iE2t )v a + cos 9 sin 9(-e~ iElt + e~ im )vp. (4) 
The probability of finding the original flavor, referred to as the survival probability, is 

Pu a ^u a = I < v a \v a {t) > | 2 = 1 - sin 2 (2#) sin 2 ( 2 ~ 1 tj 

= 1 - sin 2 (2fl)sin 2 ( 1.267 Atti 2 (eV 2 ) /^n ) ( 5 ) 

and the probability of finding the other flavor, called the appearance probability is 

P^ vp = sin 2 (2#) sin 2 (l.267A7n 2 (eV 2 )-^|^ y ) . (6) 

Here we take the approximation Ej ~ \p\ + 7nj/2E u and denote Am 2 = m 2 — m 2 . 
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The characteristic behavior of this expression as a function of sin 2 (26) and Am 2 is the following: 
for large Am 2 , the argument of the sine function is large and hence oscillates rapidly in even a very 
small energy range. The energy average of the sine function involved becomes ^, hence we have 

W(26)^2P Ua ^. (7) 

For small P Va ^ V(j and near sin 2 (26) = 1 the curve behaves like 



2 \ 



Am 2 « -i — (8) 



Hence the effective probe in Am lies in the region 



< Am < (9) 



1.267 L/E u 1.267L/E U 

This feature is true even for cases of more than two flavors of neutrinos. If the different mass 
squared differences are well separated, the effective probe of each mass squared is in the region of 
L/E u which satisfies the above relation. 



2.2 Evidence of neutrino oscillations 

The most convincing evidence of neutrino oscillations is given by the Super-Kamiokande (Super-K) 
experiment on atmospheric neutrinos ||. The solar neutrino data provide the corroborative evi- 
dence of the oscillatory scenario. The reactor experiments have not seen any evidence of oscillations, 
but they provide strong constraints on some of the allowed regions of the oscillation parameters 
and are complementary to the atmospheric and solar data. The accelerator based short baseline 
experiments, LSND and KARMEN, are so far inconclusive. Fig. |l| summarizes the pertinent results 
of all experiments to date. The filled areas are allowed sin 2 (26) — Am 2 regions that provide the 
evidences of neutrino oscillations. Excluded regions lie above and to the right of the corresponding 
curves shown. Sensitivity curves of some of the future experiments are also shown. 

2.2.1 Atmospheric neutrinos 

The atmospheric neutrino data are the so-called smoking gun of neutrino oscillations. They offer 
the strongest evidence of neutrino oscillations to date. The first indication of the "atmospheric 
neutrino anomaly" was observed in the 1960's by experiments in South Africa and India and 
was confirmed during the 1980's to 1990's by 1MB, Kamiokande, FREJUS, NUSEX, Soudan and 
MACRO experiments. The data are usually analyzed in terms of the double ratio, namely the 
observed over expected ratios of muon to electron neutrinos, to reduce the systematic error. Such 
ratio of ratios is measured to be about 0.60 ± 0.06 instead of 1 as expected, indicating neutrino 
oscillations. The breakthrough came in the 1990's when the Super-Kamiokande experiment went 
online. With a much larger statistical sample and broad energy range, they observed the following 
effects specific to neutrino oscillations ||: 
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10" 4 10" 3 0.01 0.1 1 



sin 2 20 

Figure 1: Summary of current experimental results on neutrino oscillations. Filled areas are allowed 
regions from the evidence of neutrino oscillations. Excluded regions are up and right of curves. 
Sensitivity curves of future experiments are also shown. 
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Zenith angle dependence: The down-ward(# = 0°) and upward(# = 180°) going neutrinos 
travel a distance varying from L=20 km to L= 13000 km. The observation of zenith angle de- 
pendence of muon neutrino events consistent with the expected variation of oscillation probability 
as a function of L. For electron neutrinos, such variation was not observed and it is consistent 
with reactor neutrino experiments which found null oscillations of v e |l^, [l3|]. Such zenith angle 
dependence of has also been observed later by the Soudan 2 and the MACRO experiments. 

Energy and L/E dependence: Muon neutrino events as a function of the neutrino energy 
and L /E have been observed and they are consistent with expected variation of neutrino oscillation 
probability. Again, electron neutrinos have no such variation. 

East-west anisotropy: The Earth's magnetic field, which deflects charged particles traveling 
in the atmosphere, affects the cosmic charge particles moving eastward differently from those moving 
westward. The neutrino flux model which takes into account the effect of the Earth's magnetic 
field, is consistent with observed east-west anisotropy for horizontal going //-like and e-like event, 
supporting the neutrino oscillation interpretation 

The Super-K data allows for the — * v T oscillation in a region shown in Fig. [j] at 90% CL. It 
is also reported recently |l5| that the muon neutrino to a sterile neutrino, — ► u s , as a dominant 
oscillation mode is excluded at the 99% confidence level (CL). Another exotic mechanism, i.e., 
the neutrino decay, although not completely excluded, is not a satisfactory explanation for the 
atmospheric neutrino data 0]. 



2.2.2 Solar neutrinos 

The "solar neutrino deficit" was first detected in 1968 by the Homestake experiment and 
later confirmed by the Kamiokande, GALLEX, SAGE and Super-K experiments. The observed 
neutrino flux is only 1/3-1/2 of what was expected by the Standard Solar Model (SSM) depending 
on the neutrino energy. During the last decade, the understanding of the solar structure has been 
significantly improved, facilitating the analysis and interpretation of the solar neutrino data [|l^] . 

The only plausible explanation for this deficit is electron neutrino oscillation. Fig. [l] shows 4 
allowed regions for v e oscillation into z/„ or v T in the Am 2 — sin 2 (20) space from results of all solar 
neutrino experiments: 

Vacuum long wavelength oscillation solution (VAC): This solution assumes no matter 
effect and is also known as the "just so" solution. The best fit gives a very small squared mass 
difference and nearly maximal mixing: 

Am 2 ~ 8 x 10~ n eV 2 , sin 2 (26>) ~ 0.8. (10) 

In this solution the Earth-Sun distance is approximately half of an oscillation length, L osc = 
4irE u /Am 2 , which implies a significant dependence of the v e survival probability on the energy of 
the neutrino. 

Small mixing angle MSW solution (SMA): In this solution the matter effect is included. 
The best fit of the data prior to the most recent Super-K data gives: 

Am 2 ss 5 x 10" 5 eV 2 , sin 2 (26») ps 10~ 3 - 10~ 2 . (11) 
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Large mixing angle MSW solution (LMA): Again the matter effect is included and the 
allowed parameters are: 

Am 2 ss 2 x 1(T 5 eV 2 , sin 2 (26>) ss 0.8. (12) 

Low mass solution (LOW): This is similar to the LMA case but with a much lower squared 
mass difference. The parameters are: 

Am 2 ss 8 x 1(T 8 eV 2 , sin 2 (26>) « 1. (13) 

For sterile neutrinos, only VAC and SMA solutions are allowed at 99% CL with similar param- 
eters given above. However, recent Super-K results [19] have ruled out both the sterile neutrinos 



and the VAC and SMA solutions at 95% CL. The LMA and LOW are viable with the LMA being 
the favorable solution at 90% CL. 

It should be remarked that the data do not rule out the oscillation of v e into two or more flavors 
with comparable probabilities. Furthermore, some analysis concluded that the solar neutrino data 
is still ambiguous and the SMA may still be viable pot Future data from SNO and BOREXINO 
as well as reactor data from KamLAND will be decisive for the determination of the solar neutrino 
oscillation parameters. Various possibilities of checking the solar neutrino problem can be found in 
Ref. |1. 



2.2.3 Neutrino oscillation searches at nuclear reactors 

Nuclear reactors are an intense and stable low energy source of u e , which is complementary to the 
solar neutrinos (z^ e ). The beam emitted from the reactor is isotropic with an average energy around 
3 MeV and an maximum energy of 10 MeV. The v e flux is known within 2.7% when the reactor 
parameters are given. (For references, see Ref. |22||.) 

The two reactor neutrino experiments CHOOZ |l2j and Palo Verde both have a baseline 
of about 1 km and were recently completed. Null oscillation results have been reported and they 
excluded v e — > v x oscillations in a region as shown in Fig. |l| at 90% CL. An important conclusion 
from the data of these experimental is that the atmospheric neutrino oscillation is not predominantly 
caused by the — > u e . 

The future long baseline reactor neutrino experiment, KamLAND, may find the smoking gun 
for the solar neutrino problem. The detector consists of 1 kt of liquid scintillator and can detect v e 
from all nuclear reactors in Japan with an event rate of 2/day. The characteristic baseline is 160 
km and its sensitivity curve is shown in Fig. [I]. If the solution to the solar neutrino deficit turns 
out to be LMA, KamLAND will be able to see the effect of neutrino oscillations, similar to the 
solar neutrino deficit. 



2.2.4 Short baseline accelerator neutrino oscillation experiments 

The accelerator neutrino experiments discussed below are low energy short baseline experiments 
which search for the oscillation Vuipu) — > v e {pe)- The results came from two experiments, LSND 
and KARMEN. The LSND experiment used a scintillator and Cherenkov detector with L = 29 m 
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and average neutrino energy E v = 110 MeV. It has the following results |2^] before its completion 
in 1998: 

Search for the — > v e oscillation: 

The i>n beam of LSND is produced from ji + decay at rest. A total of 22 v e events are observed 
via the reaction v e P — > e + n, while only 4.6 ±0.6 events are expected to come from the background. 
A fit to the e + energy spectrum yields an excess of 51.8j[jgg ± 8.0 events, corresponding to an 
oscillation probability of (3.ll}'J ± 0.5) x 10~ 3 if interpreted as neutrino oscillations. The allowed 
region is shown in Fig. |l]. 

Search for the — ► v e oscillation: 

In their later data sample, LSND has observed 40 events in Va — ► v e by using from n + decay 
in flight. The v e signal is detected via the charged current reaction v e C — ► e~ X. The expected 
backgrounds consist of 12.3 ±0.9 events from cosmic ray and 9.6 ± 1.9 from other neutrino induced 
processes. Interpreted as neutrino oscillation, the excess of 18.1±6.6 events correspond to a i>u — > v e 
probability of (2.6 ± 1.0) x 10~ 3 which is consistent with the i>n — ► v e oscillation probability given 
above. Since this is a 2.6<r effect, confirmation with higher statistics is necessary. 

The KARMEN experiment uses a liquid scintillator detector with L = 17 m at a comparable 
neutrino energy. KARMEN has not observed any oscillation signals. Its excluded region is also 
shown in Fig. [l|. The data taking by KARMEN will be finished in 2001. Its sensitivity to — ► v e 
will be increased by a factor 1.7, but still cannot cover the whole LSND allowed region. 

The Fermilab Mini-BooNE experiment which will be online in 2001, can cover the whole LSND 
allowed region. The Mini-BooNE Vn beam, which has a broad band of 0.3-2 GeV and v e contam- 
ination of about 0.3%, is produced by the Fermilab 8 GeV, high intensity Booster Synchrotron 
proton beam. The 445 ton mineral oil detector is located 500 m from the neutrino source. The 
expected — ► v e oscillation suggested by the LSND data will produce 1500 excess events of the 
electron type which is many cr's away from the statistical error of the experiment. The sensitivity 
curve is reported in Fig. |]. 

We summarize in Table [l] the current status on the determination of the Am 2 and sin 2 (26) 
parameters in the various types of experiments. Some of the most recent Super-K data are included. 
It should be emphasized again that the best values of the oscillation parameters are obtained 
assuming two-neutrino mixing. 

2.3 Long baseline accelerator experiments 

Below we discuss briefly four accelerator based LBL experiments: the online K2K and the three 
approved experiments, MINOS, ICARUS and OPERA. We also summarize the newly proposed 
LBL experiment J2K. 

K2K-the online LBL experiment: 

The K2K LBL neutrino oscillation experiment uses the meson neutrino beam produced at the 
KEK 12 GeV proton synchrotron. The average neutrino energy is 1.4 GeV which is below the tau 
lepton production threshold. It is expected to have 10 20 protons on target (POT) in 5 years. The 
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Exp Sources Phenomenon Recent results/remark 



Reactor 




v e survival 


No oscill observed, rule out regions 
Expected future data from: 
KamLAND 


Solar 




v e deficit 


Regions allowed: LMA & LOW & SMA & VAC 








Expected future data from: 








Super-K, SNO, BOREXINO 



Dominant (@99% CL): 

-»■ u T , v e ,v s 
Regions allowed: 

Aml tm = 3.2(1.5 - 5) x 1(T 3 eV 2 
sin 2 (2fl atm ) = 1(> 0.88) 
Expected future data from: 

Super-K, SNO 

Region allowed by LSND: 
Am£ SND = 0.2 - 1.0 eV 2 
sin 2 (26»LSND) = 0.003 - 0.03 
Controversial: some LSND allowed 
regions ruled out by KARMEN & Bugey 
Expect future data from: 
KARMEN, Mini-BooNE 



Accelerator 






survival 


Preliminary results from K2K 


hadron beam 






appear 


Three new exps online 2005 






V T ,V T 


appear 


Determine v parameters, CP(?) 


Accelerator 






survival 


New generation of exps & accelerator 


/Lt-storage 






appear 


Accurate determination of v parameters 








appear 


Search for CP/T effect 






V r ,V T 


appear 





Atmospheric v e deficit 
smoking gun 



Accelerator v^v^ u e , v e appear. 
SBL 



Table 1: Summary of neutrino oscillation experiments. 
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physics goals of the K2K are: (a) to check the v„ survival probability, (b) to ensure that there is no 
large v e appearance probability, and (c) to determine the oscillation parameters. The far detector 
is the Super-Kamiokande located 250 km west of KEK. The experiment was commenced in April 
1999 and the first neutrino event was observed on June 19, 1999, formally established the feasibility 
of this new type of high energy physics experiment. 

Results from the first year of running corresponding to 2.29 x 10 19 POT, were reported at 
ICHEP2000 |24}| . A total of 27 interactions were observed inside the central detector, while 
40.3^4g were expected from the measurement of the front detector if f„ is not subject to oscillation. 
The data is consistent with i/„ oscillation with Am 2 approximately 1 x 1CP 3 eV 2 , in agreement 
with the atmospheric result of Super-K. However the statistics is not significant enough to draw 
any firm conclusion. Nevertheless, it can be said that the data disfavor null oscillation at the 2a 
level. When more data are accumulated, an oscillation analysis of the neutrino energy spectrum 
can be performed, a stronger conclusion on oscillations can be made and the favorable parameter 
region may be obtained. 



MINOS[25|: The neutrino beam comes from the 120 GeV Fermilab proton synchrotron and 
the detector is located in Soudan, Minnesota, USA. The distance between the beam source and the 
detector is 730 km. A brief summary of the experiment is given below: 

• Goals: To (1) detect — > u T , (2) measure Am 2 tm and sin 2 (2# atm ) to 10%, and (3) search for 
the Vn — * v e component of the oscillation. 

• Detector: 5 kt Iron-scintillation sandwich calorimeter with toroidal magnetic fields in the thin 
steel plates. 

• Expected number of charge current interactions: 

Beam regime energy (GeV) CC events/kt-yr 
high 12 30000 

medium 6 1450 

low 3 450 



ICARUS pU The neutrino beam CNGS (CERN Neutrino Beam to Gran Sasso) is derived from 
the CERN 450 GeV proton synchrotron and the detector is located in the Gran Sasso underground 
laboratory. The average energy of the neutrino is less than 30 GeV and the distance traveled by 
the neutrino is 743 km. 

• Goals: To do both LBL and atmospheric neutrino oscillation experiments. It can observe all 
three channels of neutrinos: v e , and v T and is optimized for the observation of v T . 

• Detector: The detector is a liquid argon image detector. It involves new detector technology 
- liquid TPC - with modular structure. The first module of 600 t will be installed soon for 
solar neutrino and atmospheric neutrino physics. 

• Excellent e identification for the v e appearance experiment. 

• v T appearance is a key physics goal. ICARUS expects 600 v T in the liquid target. For 4 
years running the experiment sensitivity will be increased to Am 2 tm = 1.3 x 10 -3 eV 2 and 
sin 2 (26»23) = 1.2 x 10~ 2 . 
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OPERA p7| : The beam profile allows the CERN CNGS neutrino beam to supply neutrinos to 
both ICARUS and OPERA which is also located in Gran Sasso. 

• Physics Goal: Optimized for — > v T oscillation and excellent electron identification for the 

— > v e search. 

• Detector: Super module constructed out of individual scalable modules. Emulsion cloud 
chamber construction: massive dense material (Pb/Fe) plates as target plus thin emulsion 
sheet to track the r decay, target section followed by a muon detector to reduce background 
to a very low level. 

HIPA to Super-K (J2K) J28|: A new proposal has been made to use the proposed nuclear 
physics facility, the High Intensity Proton Accelerator (HIPA, formerly the JHF), to generate 
an intense neutrino beam with Super-K as the detector. HIPA, located about 60 km north-east 
of Tokyo, is a high intensity 50 GeV proton synchrotron accelerator, with an intensity of 1 x 
10 21 POT/year. It will take 6 yeas to complete starting from 2001. We summarize the physics goals 
of J2K in the following: 

• Observe Am 2 to an accuracy of 2 x 10 -4 eV 2 . 

• Assuming the dominant oscillation — > u T , sin 2 (28^) = cos 4 (#13) sin 2 (2023) will be measured 
with an accuracy of 0.04. 

• v e appearance search: to measure sin 2 (2# e ) = sin 2 (#23) sin 2 {28 1^) to 0.05. 

• v T appearance search: Very little hope unless Am§ 2 is much larger than the current value of 
3.5 x 10" 3 eV~ 3 . 

• search for the possible presence of sterile neutrino v s by neutral current events can provide a 
stringent limit on the existence of v s . 

Note that because of the distance, the matter effect is not significant and the oscillation effect is 
similar to that of the vacuum. 



15 



3 Theoretical introduction to neutrino oscillation 



The experimental data available to date have provided the information needed for the construction 
of a generic framework for massive neutrinos. Let us first summarize the relevant experimental 
data before we embark on describing the possible theoretical scenarios that embody the data. If 
we accept all data as discussed in the proceeding section we see that there are three distinctive 
mass scales provided by the three categories of experiments: the LSND, atmospheric, and solar. 
The mass square difference (MSD) and the mixing angle in each category of experiments are given 
in Table We list the current best values which are mostly from the Super-K collaboration 
for the atmospheric and solar neutrinos and the LSND collaboration for the SBL experiments. 
Coming from a single experiment, the LSND data have to be considered as tentative and require 
confirmation. 



Category of Exp 


MSD (|Am 2 | (eV 2 )) 


mixing angle (sin 2 (29)) 


LSND 


0.2-1 


0.003-0.03 


Atmospheric 


(1.5 - 5) x 10~ 3 


> 0.88 


Solar LMA 


2 x 10~ 5 


ps 0.8 


SMA 


5 x 10~ s 


0.001-0.01 


LOW 


8 x 10~ 8 


PS 1 


VAC 


8 x 10 _il 


ps 0.8 



Table 2: Summary of mixing parameters of the three categories of neutrino oscillation experiments 



The striking feature of the MSD pattern given in table [2] is that the three category of experiments 
provide three well-separated mass scales. This hierarchical structure of neutrino masses is similar 
to that of the quarks. It requires four distinct masses and hence the existence of at least four 
neutrino flavors. Therefore, accepting the LSND data immediately implies a non-trivial extension 
of the neutrino sector of the SM. The additional neutrino not contained in the SM is devoid of 
interactions with the SM electroweak gauge bosons, usually referred to as the sterile neutrino and 
denoted as v s . If the LSND data are excluded, the three SM neutrino flavors are sufficient and 
no extension of the number of neutrinos is necessary. The latter case gives rise to the 3-neutrino 
scenario and the former the 4-neutrino scenario. Within either scenarios there are several cases 
which are different from one another by the ordering of the masses of the individual neutrino mass 
eigenstates. We will illustrate the different mass assignments below. 

In view of the uncertainty of the LSND data, the discussion will be focused on the 3-flavor 
scenario. Comments on the 4-neutrino scenario will be made at the end of this section. For a 
detailed introductory theoretical review we refer to Ref. 2S]. 

We define the neutrino states and masses by: 

Flavor eigenstates: v a , up, v^... 
Mass eigenstates: ui, U2, v-$... 
Masses: mi, m^, m^... 
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The column vector of the flavor states will be denoted as [i> a ], and the column of mass eigenstates 
by [uj]. The two sets of states are related by a unitary matrix U: 

["a] = U[vj] 
U = (U aj ) (14) 

To illustrate the different cases of the mass assignment let us consider the 3-flavor scenario. 
Because of the order of magnitude difference in their MSD, the solar and atmospheric data imply 
that two of the mass eigenstates lie closely in their mass values which we will assumed to be V\ 
and 1/2, with the third, assumed to be ^3, relatively far separated from the first two. Since the 
existing data give only the values of the MSD's not their signs, the mass eigenstates cannot be 
unambiguously identified. Let us denote the mass order mj < m| < mf as the jkl-case, then 
there are four possible cases of the mass orders: 123, 213, 312 and 321. These four possibilities 
corresponding to four possible sign assignment to Am^ = m\ — m\ and Am^ = m| — m§. 
Future experiments will have to find out which case is correct. Three of the four possible level 
assignments, 123, 312 and 321 are shown in Fig. [3| 



2 1 

j A?7ig 0lar Amg 0lar 



Am atm 



^ Amg 0lar 



Am atm Am atm 



Figure 2: Three possible level structures of the 3-neutrino scenario 



3.1 Oscillation in vacuum 

The formulae given in this subsection are valid for an arbitrary number of neutrino states. 
3.1.1 Oscillation probabilities 

The oscillation probabilities are functions of the mixing matrix elements U a j, the MSD's Am| - = 
m| — m| with m\ > m|, the neutrino energy E u , and the oscillation length L. The oscillation 
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probability of u a — > up is given by 



Pu a ^u fi = S a(3 -Aj2Y^sm 2 (A kl )-2j2J 3 a^m(2A kj ) (15) 

j<k j<k 



where 



Y* ee Re(U aj U* k U* pj Up k 
J J i = Im(u aj U* ak U* pj U pk ) (16) 
A kj = 1.267A<. ( eV 2 )^. 

Y^p is symmetric in j/c and a(5, is anti-symmetric in jk and a/3: 

— v k i — v3 k — V k 3 
1 a/3 — 1 f3a ~ 1 /3a ~ 1 a/3 

Ja(3 = ^ = -41 = -^ ( 17 ) 

CP and T violation effects will be exhibited if J^L 7^ 0. 

The corresponding anti-neutrino oscillation probability can be obtained from the above expres- 
sion by the replacement U — > U* , giving 

Pv^vp = <5^-4^y^sin 2 (A fcj ) + 2^jSsin(2A fc ,) (18) 

j<k j<k 

This replacement for the anti-neutrinos together with the symmetry relations of Y^ and given 
above allows us to write down all the probability formulae for a given pair of neutrinos, u a and up. 

3.1.2 CP/T and CPT asymmetries 

For a given pair of neutrino flavors, the probability P Va ^ Vp is related to three other oscillations 
probabilities by CP, T and CPT transformations: 



CP : 


Py a 


^rVp — 


* Pfa 




T : 


P 


-rVf} — 


> P 




CPT : 


P 


^rVf} — 


► P- 





(19) 



Effective measurements of symmetry violations can be made by the so-called asymmetries as 
defined below. For a given pair of flavors of neutrinos and anti-neutrinos, v a , up, D a and Dp, six 
asymmetries can be defined for u a ^= up: 

CP asymmetries: 

>cp _ Pyg^yy ~ Pyoc^vp _ ^2j<k jjp sin(2Afcj) 

A a/3 - p , p_ _ - jk . 2 . . ) l zu J 

^v a -^up + ii/ a -i/ fl 2 Y, j<k Y J a p sm z (A ifc ) 
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T asymmetries 

- p —5 W 

CPT asymmetry 

A «/3 - p T~5 

Three more asymmetries, one each corresponding to the above asymmetries, are obtained by 
the interchange of v a and up. 

The analytic expression of the oscillation probability assumes necessarily CPT symmetry and 
therefore gives vanishing CPT asymmetries. Hence given CPT symmetry there is only one inde- 
pendent asymmetry because of the identity, 

A%p = -Ap£ = A T afi = -Ap a with CPT invariance in vacuum (23) 

due to the symmetry properties of Y^p and J^. However, the CPT asymmetries given above offer 
a simple way to check the underline assumptions of the CPT invariance in the neutrino sector and 
should be made. 

In a neutrino factory, all six asymmetries can be measured for the and v e channels. But 
in the case of conventional meson-neutrino beams only the CP-asymmetries A^ and A^f are 
accessible. 



3.2 The three-flavor scenario in vacuum 

The three-flavor scenario is the experimentally favorable scenario. Although matter effects are 
always present in LBL experiments and their inclusion is necessary in order to extract precisely 
the oscillation parameters, the consideration of the vacuum case is nevertheless useful. It is a 
good approximation for experiments with shorter baselines. Moreover, it gives relatively simple 
expressions for the various oscillation probabilities which provides a more transparent picture of 
the physics. This allows us to look for the optimal experimental conditions to conduct a particular 
LBL experiment. 



3.2.1 Basic formulae 

The unitary mixing matrix can be parameterized as: 



U 



I 

V 



C12C13 
-C23S12 - C12S13S23 
S12S23 - C12C23S13 



C13S12 
C12C23 - S12S13S23 
-C12S23 - C23S12S13 



'13 



where Sjk 



sin(# 



jk)i Cjk 



COS 



tjk) 



and s 



j fc 



sin 



fcje 



C13S23 (24) 
C13C23 / 

l5 . 6jk defined for j < k is the mixing angle 



of mass eigenstates Vj and and 8 is the CP phase angle. The presence of a non-vanishing S will 
give rise to CP- and T-violation phenomena. 
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The CP/T- violation effect is given by the Jarlskog invariant |3(j]. From the unitarity condition 
it can be shown that there is only one Jarlskog invariant in the three-neutrino scenario. Hence 
the CP/T-violation part of all the oscillation probabilities are the same except for a possible sign 
difference. The Jarlskog invariant, denoted as J (5) is defined by 

J af} = J ( 6 )J2 € ^7 e jkl ( 25 ) 

1,1 

where (a/37) i n the e-symbol run over the three neutrino flavors with 6 e ^ LT = +1 and other permu- 
tations conventionally defined, and 

J(6) = -ci3Sin(20 12 )sin(20 13 )sin(2^3)sin(J). (26) 
o 

The CP-violation term can be rewritten as 

" 2 E J S sin ( 2A fci) = " 8J (<5) sin(A 21 ) sin(A 32 ) sin(A 3 i) £ e Q/37 (27) 

j<k 7 

The explicit oscillation probabilities can be written done straightforwardly, including those for 
anti- neutrinos with the replacement U a j — ► U*j. 

3.2.2 Identification of mixing angles and MSD 

In terms of the experimental data, the mixing and MSD parameters can be identified as follows: 

Am| 2 = Am 2 atm 

#12 = Qsolar (28) 
#23 = #aim • 

The ranges of values of the above parameters are given in Table || summarizing the experimental 
data given at the beginning of this section. The MSD Am^ is not independent, 

AmJi = Am§ 2 + Am 21 . (29) 

The third mixing angle, #13 is not precisely measured. The CHOOZ collaboration gives an upper 
bound of the angle, 

sin 2 (26»i 3 ) < 0.1. (30) 

3.2.3 Mass hierarchies, regimes of low and high mass scales 

From the hierarchical structure of the MSD 

lAmlJ « |Aro§a| » |Am^|, (31) 
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we can classify LBL's according to the neutrino energy and the baseline length. Let us define the 
MSD scale as the smallest MSD that gives the maximal oscillation, i.e., Ajk ~ vr/2 

A-Lc(eV 2 ) - l-240^p (32) 

4c ~ and A ™osc ~ 

Am?, 



which suggests that there are two relevant experimental regions: Am„ sc « Am^ and Am^ 



Low mass scale regime Am^ c rs Am^ 

This region corresponds to the solar neutrino oscillation regime and for low E u and long L, 
where sin 2 (Am 3 i) and sin 2 (A32) oscillate rapidly and can be replaced by \. Then the neutrino and 
anti-neutrino oscillation probabilities can be simplified as 

P Ua ^ = 6 af3 (l - 2|*7 a3 | 2 ) + 2|[/ a3 | 2 |£/ /33 | 2 " * Y a$ sin 2 (A 21 ) - 2J* 2 sin(2A 21 ) 
P Pa ^ = S af3 (l - 2|*7 q3 | 2 ) + 2|f/ a3 | 2 |^ 3 | 2 - 4Y$ sin 2 (A 21 ) + 2J* 2 sin(2A 21 ) (33) 

In this mass scale regime, unless the baseline is very long, the neutrino energy will have to be 
low. The experimentally interesting measurements are the v e survival probability as in the solar 
neutrino oscillation experiment, and the — > v e appearance probability: 

Pv^v e = l~\ sin 2 (20 13 ) - 4 sin 2 (2fl 12 ) sin 2 (A 21 ), 

P v ^ Ve = ^ 2 3 sin 2 (2^ 13 ) + c 2 3 4sin 2 (20 12 )+2J(<5)sin(2A 21 ). (34) 

The v e survive probability measures sin 2 (2#i 3 ), Am 21 and c 4 3 sin 2 (2#i 2 ) if the ratio L/E v can 
be varied so that the constant term and the oscillating term can be separately measured. 

In this vacuum approximation the favorable distance to neutrino energy ratio for the LMA 
solution of the solar data of Am 2 olar = 2 x 1(T 5 eV 2 is L/E v = 6.3 x 10 4 km/GeV. For a neutrino 
energy of 3 MeV, which is in the average energy of the neutrino generated in a reactor, the distance 
is 190 km. As discussed in the preceding section, this is in the range of the KamLAND reactor 
experiment. 

Note that the CP-violation term in the — > v e oscillation is generally small as it is proportional 
to sin(20i 3 ). 

High mass scale regime Am 2 sc « Am 31 « Am 32 

This region is suitable for the terrestrial and accelerator neutrino oscillation experiments. In 
the leading approximation, the oscillation probabilities can be approximated by dropping the terms 
proportional to Am 21 . We will also ignore the CP-violation term in this approximation but will 
consider at the end of this discussion. 

~ ^-4|[/ a3 | 2 (l-|f/^ 3 | 2 )sin 2 (A 32 ). (35) 
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In more detail we have 











a 1 - 4c 2 3 s 2 , 3 (l - c? 3 s| 3 ) sin 2 (A 32 




-*V e ? 






s l-sin 2 (2^3)sm 2 (A 32 ), 






- *v 


-»J> e p 


a s 2 23 sin 2 (2%)sin 2 (A 32 ), 








-+Z> T ? 


»4 sin 2 (20 23 ) sin 2 (A 32 ), 


Pu e - 






■+£y p 


^4 sin 2 (2%)sin 2 (A 32 ). 



(36) 



Due to the smallness of sin 2 (2^13), in this experimental regime the appearance probabilities 
P u -+ Ve and P Ve -^ Vr are small, probably at a few percent level. However, because the mixing angle 
# 2 3 is near maximal, the appearance probability P v -* Ut can vary from close to 1 to zero when 
L/E v varies and it can provide a good measurement for the product cf 3 sin 2 (2# 2 3). The survival 
probability for v e stays large. But the survival probability for varies from close to 1 to near zero 
when L/E u varies. 

Similar to the case of the low mass scale regime, the CP violation term is proportional to 
sin(2A 2 i) and is neglected in the above approximate expressions. But it can be identified as: 



AP° — P — P- 



->CP 

» +4/(5) sin(2A 21 ) sin" (A 31 ) (37) 

and the corresponding CP asymmetry, 

A CP „ sin(2g 12 )ci3C 23 sin(2A 21 ) 
^ 2s 2 3 sin(20i 3 ) 

This indicates that in order to be able to see the CP-violation in the Vn — > u e oscillation in this 
high mass scale regime, the solar mass scale Am 21 can not be too small and it is favorable to have 
a large L/E v to maximize the product sin(2A 2 i) sin 2 (A3i). 

Identifying Am 2 sc ~ |Am 2 2 | = 3.2 x 10 _3 eV 2 , then the effective experimental probe for a 
baseline length of 2100 km gives E v ~ 5.4 GeV. This shows that the neutrino energy between 1 
to 10 GeV from the H2B neutrino beam is in the optimal range for the probe of the atmospheric 
oscillation mass scale. 



3.3 Oscillation in matter 

The vacuum oscillation formulae is modified by the presence of matter along the path of the 
neutrino. When neutrinos propagate through matter, the SM neutrinos can interact with the 
quarks in the nucleons and electrons in the atoms to undergo both elastic and inelastic scatterings. 
The inelastic scattering and the elastic scatter off the forward direction will cause attenuation of 
the neutrino beam. Since the cross sections are extreme small, the attenuation are insignificant. 
However, the elastic scattering in the forward direction is a different matter. Although it does not 
change the direction of the neutrinos in the beam, it can modify the vacuum mixing angles and 



mass eigenvalues of the SM neutrinos. This matter effect is the well-known MSW effect [31], which 
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is not CP symmetric, hence the modifications to the anti-neutrinos are different from those of the 
neutrinos. The matter effect on the electron neutrino is different from that on the muon or the 
tau neutrinos. It should also be noted that the matter effect is cumulative, in analogy with the 
propagation of light in a medium; the effect will be manifested more clearly when the length of 
propagation increases. Hence a clear detection of the matter effect will require a sufficiently long 
baseline. This intuitive conclusion is born out by the explicit calculation on the effect of the matter, 
and the H2B's 2100 km baseline can provide explicit checks of the matter effect. 

For the electron neutrino scattering through the matter, there are both charge and neutral 
current interactions. The muon and tau neutrinos subject only to neutral current interactions. For 
a sterile neutrino, since it does not interact with the SM gauge bosons, its propagation will not be 
affected by the presence of the matter. There are special considerations to simplify the calculation 
of the matter effect. We give some more details below. In the case of three-neutrino scenario 
without a sterile, the neutral current interaction is the same for all flavors. The effect of the elastic 
forward scattering due to neutral current is to give a common phase to all flavors. This has no 
effect on the oscillation and the neutrino current interaction can be ignored. However, the charge 
current interaction, which contributes only to the elastic scattering of electron neutrinos with the 
atomic electrons, has to be considered. Because the ev e and ev e scatterings are different, matter 
effects on neutrinos and anti-neutrinos are different. Hence, the matter effect has to be carefully 
removed before the information on CP-violation can be extracted. The effect of T-violation can be 
investigated more readily, independent of the matter effect. 

In a LBL experiment, as the neutrino traverses through the Earth, it encounters the Earth 
matter which may vary along the neutrino path. The density dependent matter effect is commonly 
dealt by the Schrodinger equation approach. The time evolution of the neutrino is equivalent to, 
and can be expressed as a differential equation in the distance of its propagation, 



.d_ 

'dL Ua ~ 2K 



\ £ ((£ ™PajUh) + ASaetp^j vp (39) 



A = 2E u a c 

a c = V2G F n e = V2G F N A Y e p 

where ac is due to the charge current effect, Gf is the Fermi constant, n e the number density of 
the electron, Na the Avogadro's number, p(gm/cm 3 ) the matter density in units of gram per cm 3 , 
and Y e the average of number of electrons per nucleon. The matter density n e and hence p depend 
on L. For most of the Earth density, Y e can be taken as 1/2. 

The above equation can be integrated numerically for the solution. Note that as mentioned 
before we have dropped the neutral current effect, which is common to all the three SM neutrinos 
and can be ignored. The expressions for anti- neutrinos can be obtained by the replacement: U —* U* 
and A — > -A. 

For constant matter density, analytic expressions of the mixing parameters and oscillation 
probabilities expressed in terms of those of the vacuum quantities exist for the three-neutrino 
scenario ]32|, [33|] . The explicit expressions for constant matter density are useful in estimating the 
magnitude of the matter and CP/T effects. Especially the approximate expressions are transparent 
in their physical meaning. The rest of this subsection is devoted to analytic expressions for the 
case of constant matter density. 
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3.3.1 Approximate expressions for constant matter density for high mass scale 

In the following we list some relevant approximate expressions in the case of constant matter density 
in the leading oscillation for Am^ c ~ Am^ ~ Am^ >> Am^, which are relevant to terrestrial 
neutrino oscillation experiments and those for L/E u of the order of (Ami) -1 , such as the H2B. 
Quantities that are subjected to modification by the approximated matter effect will be marked by 
the superscript "m" . The approximate expressions are simple enough that their physical meaning 
is clear. 

A comparison of the approximate expressions with the exact one will be made at the end of the 
next subsection. The expressions given below, that set both Am^ and the CP phase 5 to zero, are 
taken from Ref. f|327| 



pm 

f- L ~ 




« 4sm 2 (2^3)sm 2 (A^ 2 ) 


1 

~ 4 


pm 




« sin 2 (20 23 ) (sin 2 ^) sin 2 (A£) + cos 2 ^) sin 2 (A£) 


pm 




_ i pm pm 




pm 

Ue~ 


->V T 


« c 2 3 sin 2 (2^3)sin 2 (A™) 




pm 




= l- S in 2 (2^3) S m 2 (Ar 2 ) 




where 









(40) 



-*»<»5) = 2^ (41) 



with 



S = V(cos(2^i 3 ) - A/Am^Y + sm 2 (29 13 ) (42) 
A^SAo, W = \(l + t±r + 8)Ao, A^\{i + ^- S )a^ (43) 

A » s ^ = ° 67A " i - (eV 'liciV) (44) 

The matter effect can be written explicitly for use in numerical simulation. A brief discussion 
of the earth density profile can be found in Sec. 4.1.3. 

A = 0.7634 x 10" 4 p(g/cm 3 )y eJ B 1 ,(GeV) (45) 

where Y e = 0.5 is the number of electron per nucleon. Note that the matter resonance enhancement 
occurs at 

E» « 15(GeV) ( - rrS" ' ) [ ) cos(20 13 ) (46) 



,, (r I Amj^eV 2 ) \ f l.5(g/cm 3 ) ^ 
15(GeV) x 10-3(eV2) ) [ P (g/c^)Y e ) C0S{29 



For the H2B experiment, taking p ~ 3 gm/cm 3 , then the matter resonance occurs at E v ~ 
15 GeV which is in the range of the HIPA neutrino beam. 

As already stated that the expressions of the corresponding anti-neutrino oscillation is obtained 
by the replacement A — ► —A. 
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Because of the smallness of sin 2 (26is) the above probability expressions show that the Ua — ► v e 
and v e — ► v T appearance probabilities are small, similar to the value in the vacuum case as discussed 
early. But they are enhanced due to the matter resonance effect. The z/„ — > v T appearance 
probability as well as the and v e survival probabilities are large. 



3.3.2 Exact results for constant matter density 



In the following we list the exact expressions of 3-flavor mixing for constant matter density [33|. 
Quantities modified by the exact matter effect will be denoted by a superscript or subscript "(m)". 
We rewrite the Schrodinger equation as 



jr ( m )i 
' dL 
Hm 



H v + 



(m)i 



2E„ 



-AL 



-Lu(AM 2 + AU^I ee U)W 




(47) 
(48) 



The matter modified mass matrix, denoted by O can be diagonalized by a matrix V, 



n 



AM 2 + AU ] I ee U f 
VMf m) V^ 



where 



Mi 



in) 



( Ai 



V 



(49) 
(50) 



A? 



contains the new mass square eigenvalues with the mass square eigenvalues [ 33 j 

1 



Ai 
A 2 
A3 



x 2 - 3y [z + J 3(1 - z 2 ) 



x 2 -3y[z- J3(l - z 2 ) 



(51) 



x + 3y x 2 — 3y 



where 



x 

y 



Amli + Amfi + A 

Am^Amli + A [Am 2 21 (l - \U e2 \ 2 ) + Am 2 31 (l - \U e3 \ 2 ) 

(l 1 f 2x 3 - 9xy + 27AmhAmt\U el \ 2 \\ 
cos - cos " 11 



(52) 



2(x 2 - 3yf/ 2 
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The elements of the diagonalization matrix V are given by [p3| 

A 



V » Dj ' 



V jk = Trk (\ j -Amf 1 )U: j U ek 



(53) 



where j ^ k / I takes the values 1, 2 and 3 and 



D 



in 



Hi 



N 2 + A 2 \U ej \ 2 



A 



(A j - Am 2 kl ) 2 \U el \' z + (Xj - Ami x f\U t 



(Xj - Am 2 kl )\U el \ 2 + {Xj - Aml)\U ek \ 



2 \2i 



ek 



(54) 



Now the Hamiltonian with the matter effect can be rewritten in the same form as the vacuum 
case: 



H 



M 



J^U^M 2 m) U^ 



uim) = {u £) ) = uv= ^ UakV \ 

TT im) _ (-I)'- 1 



Da 



NjU aj + AU ej (Xj - AmiMUd + (Xj - AmfJU^U, 



(55) 
(56) 



Now the matter case can be written in the same expression as the vacuum case with the matter 
modified mixing matrix and the corresponding MSD's, 



A (m) m^ = -yV- 3^3(1-^) 



A^m 2 31 



y^-ty (a* - ^3(1 - 



(57) 



-^ 2 -32/(3z + ^/3(l-z 2 ) 



We define the oscillation argument for j < k: 



(58) 



For the anti-neutrino, the corresponding quantities are obtained as in the neutrino case by the 
replacement: 



Tji m ) 

aj 
^kj 



U. 



(m) 
aj 

A ( m ) 
^kj 



uff*(A^-A), 
^kf( A ^- A ) 



(59) 
(60) 



The oscillation probability expressions are similar to the vacuum case, 



P, 



(m) 



Sap ~ 4 £ ** (uifu^U^U^) sin 2 (A 
j<k 



kj i 
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-8J^(5)J2^Y[ s[ ^ A ff) (61) 

7 j<k 

j<k 

+8J^(5)J2^U s[ ^ k 7 } ) (62) 

7 j<k 



where 



j^(5) = L A fl A ™h A Th 2 W) (63) 

j(™)(<5) = j("0(,J)(A-»-A) (64) 

J(5) is the vacuum expression given before. 

Unlike the vacuum case, the difference P^\ Vfi - Pp^Lug contains both the CP-violation ef- 
fect and the matter effect. However to estimate the CP angle we can calculate the T- violation 
asymmetry in which the matter effect and the CP phase can be separately isolated, 

P^l V(i - P^\ Va = -16j( m )(5)sin(A^ ) )sin(A^ ) )sin(A^ ) )^e^ 7 , (65) 

7 

In the approximation of the last subsection, the T-symmetry is then given by 

A (m)T 



p(m) p (m) 

pW I p( m ) 



%J (m \b) sin(A^7)sin(A^ } ) 
S | 3 sin 2 (2^3)sin(A( m) > 



(66) 



»32 J 

which can be measured at a neutrino factory but not with a meson-neutrino beam. 

Let us comment briefly on the validity of the approximate expressions given in the preceding 
subsection. In the H2B region, i.e., L = 2100 km and E v in the range of 1-10 GeV, the approximate 
expressions are good to within a few percent for the various oscillation probabilities. In the energy 
regime of a few hundred MeV and lower, the term proportional to sin 2 (A2i) is no longer negligible 
and it becomes eventually the dominant contribution. Then the approximate expressions are no 
longer valid. 



3.4 Comment on the four-neutrino scheme 

Although the current results disfavor strongly the dominance of the oscillation of v e — > u s as the 
mechanism for the solar neutrino deficit and — > u s for the atmospheric neutrino anomaly, a 
sizable contribution of the sterile neutrino to these oscillations is allowed |^ |. For a summary of 
the 4-neutrino fit of the various data, see Ref. . Because of the far reaching implications of such 
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a scenario, it is worthwhile to maintain an interest in the scenario. If the LSND result is accepted, 
there are clear three different MSD's squared mass differences of the order of magnitude: 



A™* < eV^ 
Aml tm « 10" 3 eV 2 
Am 2 LSND « leV 2 



which calls for at least four different neutrino masses. This has the implication that a right-handed 
or sterile neutrino exists. The existing data favors the so-called 2+2 scheme. However, another 



possibility, i.e., the so-called 3+1 scheme is not completely ruled out 34]. In the 2+2 scheme the 
4 neutrino mass eigenstate states are divided into 2 groups each containing 2 levels. Within each 
group the mass separation is relatively small in comparison with separation between the groups. 
Therefore, one group gives the solar energy scale and the other the atmospheric scale. The scale 
between the groups gives the LSND oscillation. In the 3+1 scheme, the 4 states also divided into 
two groups. One group contains 3 mass eigenstates and the other group one state. The 3-eigenstate 
group provides the solar and atmospheric oscillation scales similar to the 3-flavor scenario. The 
MSD between the groups is the LSND scale. Let us label the mass eigenstates as uq, v\, V2 and z/3. 
Three possible level structures for the 2+2 scheme are given in Fig. 
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Figure 3: Three possible level structures of the 4-neutrino scenario 



The three squared mass difference are identified as: 

Am solar = |Amoi| 
A ™atm = l A ^32l 

Am LSND = l Am 3ll ~ l Am 2ll> etc - 

Again the different level structures have different signs for Am?,-, j < k. In a complete deter- 



mination of the neutrino parameters, the signs of Am| • have to be determined. 
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The four-neutrino scheme is generally treated numerically with the Schrodinger equation. In 
this case there are six mixing angles and six CP phases. Three of the CP phases can give rise to 
measurable CP-violation effects in oscillation experiments. For an explicit parameterization of the 
mixing matrix we refer to Ref . || . The matter effect is more complicated than in the three-neutrino 
scenario because neutral current interactions have to be included since they are no longer common 
to all neutrino flavors. The CP effect can be sizable and may be easier to detect. 
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4 Fundamentals of LBL experiments and physics of H2B 



The neutrino oscillation is a system with a limited number of degrees of freedom, it yet exhibits a 
multitude of interesting phenomena. In the 3-flavor scenario, the system consists of 2 MSD, three 
mixing angles and one CP phase that determine the different survival and appearance probabilities. 
The numerous neutrino experiments, solar, atmospheric, reactor, and short baseline, are mostly 
looking at the survival probabilities and often use neutrino beams that exist in nature. In most 
cases there is no way to tune the neutrino beams for more desirable experimental results. Hence 
it is difficult to obtain all of the oscillation parameters for the entire mixing matrix, either because 
the statistics are too low, or the energy is not suitable. In long baseline experiments, the neutrino 
beams are produced in the accelerator according to certain physics criteria so that the experiment 
can be conducted in a controlled fashion. Ideally, the distance between the neutrino source and the 
detector can be chosen to maximize the physics output. The distance can be hundreds or thousand 
kilometers to allow high energy neutrino beams to be used and still offer suitable L/E u values. 
The LBL experiments promise to allow for detailed analysis of the oscillation parameters so as to 
provide a complete picture of the neutrino oscillations. 

We will first briefly summarize some of the fundamentals of the LBL experiment and then the 
possible physics that can be investigated at H2B. More pertinent simulation studies will be carried 
out with a specific detector later. 

4.1 Fundamentals of LBL experiments 

Long baseline neutrino oscillation experiments are not conventional high energy physics experi- 
ments. Due to the extremely weak interaction cross section of neutrinos and the long baseline 
between the neutrino source and the detector, both the neutrino beam intensity and the detector 
mass must be maximized in order to have the desired statistics. New technologies for both acceler- 
ators and detectors may be required. In the following we discuss briefly some of the fundamentals 
of the LBL experiment before we start to discuss the specific neutrino beam that may be available 
and the detectors that are necessary to achieve the physics goals. 

4.1.1 Neutrino beams 

There are two kinds of accelerator neutrino beams: the so-called neutrino factory from muon decays 
and the conventional neutrino beam from meson decays. While meson- neutrino beams have been 
built by many laboratories and the remaining technological challenge is to increase the total power 
of the primary proton beam, the neutrino factory is a completely new concept and there are still a 
host of technical issues to be worked out. 

Neutrino factory: The neutrino factory delivers a neutrino beam which contains comparable 
amount of (v^) and 9 e (u e ) obtained from the decay in a /U-storage ring. 

^~{^ + ) v^Vfj,) + v e {y e ) + e~{e + ) (67) 
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Note that the presence of both the muon and electron neutrinos is not a problem because they 
have opposite chiralities and therefore will not be a background to each other if the detector can 
distinguish between positive and negative electric charges. A description of the neutrino factory 



can be found in [36] 



The neutrino flux at baseline L from a neutrino factory of unpolarized muon of energy is 
given by 

d$ ( 2*2(3 -2x f )^ for 



dE u y 12x}(l-x f )^ for v e 

where Xf = E u /E^, uq is the number of useful decaying muon, and 7 = E^/m^ with being the 
mass of the muon. The average neutrino energies are given by 

(E Vlt ) = 0.7£ M 

(E Ue ) = 0.6£ M (69) 

Two scenarios of the number of neutrinos in a beam have been considered: no = 6 x 10 19 /year for 
an entry level factory and no = 6 x 10 20 /year for a high performance factory. For a discussion of 
the neutrino beam spread in a neutrino factor together with some sample plots, see Ref. p7j] 



Meson-neutrino beam: Neutrinos from a meson source are obtained from decays of pions 
and kaons produced by collisions of primary protons with nuclear target. The secondary meson 
beam produced by the collision is sign selected and then focused by a magnetic field. Mesons 
are then transported to a vacuum decay pipe, whose length depends on the desired energy of the 
neutrino, and finally striking a hadron absorber further downstream. The primary neutrino or anti- 
neutrino beam consists mostly of the muon flavor from ir^ and decays in the decay pipe but 
some impurities of electron neutrinos are expected due to a finite branching ratio of tt^, and 
^ decaying into electrons and positrons plus their associated neutrinos. For example, the NuMI 
muon neutrino beam at Fermilab contains 0.6% electron neutrinos. As for the energy spectrum of 
the beam, it can be a wide-band beam covering a broad range of energies, or narrow band beam 
with a selected, well-defined energy range. 

The neutrino flux at the detector site is determined by the baseline L, the number of primary 
protons on target (POT), the proton energy E p and the neutrino energy E v . The following empirical 



formula [38] describes the meson production from a proton beam on a nuclear target: 

da „ f „ d? 



X M 



2vr / x M E p -^P t dP t 
dxM J dp 

r . 1 4. K P — Dxm 

= (l-iM)' 4 (l + 5e- Dl '") (70) 



where E p is the proton beam energy, p the proton 3-momentum, xm is the Feynman x-variable 
defined as the momentum of the secondary meson divided by the momentum of proton. A, B, C, 
and D are numerical parameters which are different for different secondary particles. Table gives 



their fitted values for ir + ,ir , A' + , and K , taken from Ref. [38]. For a wide band beam, we can 
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take Eq.(^), i.e., the E^/L 2 behavior, to account for the transverse momentum spread due to pion 
decays. The neutrino flux can then be written as 

♦ ftlI ) ^"-^ + i '- C » (71) 

where we take xm — x u = 2E U /E p . It is interesting to note that this simple formula can account 
for the beam design of MINOS at various energies as discussed in Sec. 2.3. We should remark 



that there is a more complete treatment [39] for the energy spectrum for the meson- neutrino beam. 
However the difference with the above express for E v > 1 GeV is very small. Owing to its simpler 
form, we continue to use the above expression in the following calculations. 





A 


B 


C 


D 


7T+ 


2.4769 


5.6817E-2 


0.57840 


3.0894 




3.5648 


5.0673E-2 


0.68725 


5.0359 


K+ 


1.7573 


6.3674E-3 


0.81771 


5.6915 


K~ 


5.4924 


4.1712E-3 


0.89038 


2.2524 



Table 3: Numerical parameters for meson-neutrino energy spectrum 



4.1.2 Dip angle 

To deliver the neutrino beam to the detector that is L (direct) distance away, the beam has to 
point downward at a dip angle 9^ p from the horizontal plane given by 

(72) 

where te = 6371 km is the average earth radius. This is the same dip angle of the hadron decay 
pipe in the case of meson-neutrino beam. The dip angle is 1.1° for K2K, 3.3° For MINOS, ICARUS 
and OPERA. H2B has an oscillation length of 2143 km and the dip angle is 9.6°. 

4.1.3 The Earth matter density profile 

The solid earth is made of three major parts, the crust, mantle and core. The density increases 
with increasing depth D from the Earth surface. However, there are local variations of the earth 
density. A widely used model for the earth is the Preliminary Reference Earth Model (PREM) 
provided by Dziewonski and Anderson |4(| . We list in Table |4] the earth densities versus the radius 
and the depth from the Earth surface, where x = v/te, ris the radius at a given depth. The density 
is in units of g/cm 3 and the radius and depth in km. The densities given are of course the average 
densities. The extrapolation formulas are valid in the sublayers where the mass densities vary. The 
actual densities at a given depth is given in the square bracket. The earth density profile HlJ is 
also plotted as a function of the distance along a diameter from one end on the earth surface to 
the opposite end as given by the solid curve in Fig. |i] |42|]. 



,dip = Sm " fe) 
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region 


radius r(km) 


depth D(km) 


extrapolation [density] (g/cm 3 ) 


ocean 


6371 





[1.0200] ([2.6000]) 


(continent) 


6368 


3 


[1.0200] ([2.6000]) 


crust 


6368 


3 


[2.6000] 




6356 


15 


[2.6000] 




6356 


15 


[2.9000] 




6346.6 


24.4 


[2.9000] 


LID 


6346.6 


24.4 


+0.6914x [3.808] 




6291 


80 


2.6910 [3.3747] 


low velocity zone 


6291 


80 


+0.6914x [3.3747] 




6151 


220 


2.6910 [3.3595] 


transition zone 


6151 


220 


-3.8045x [3.4358] 




5971 


400 


7.1089 [3.5433] 




5971 


400 


-8.0298x [3.7238] 




5771 


600 


11.2494 [3.9758] 




5771 


600 


-1.4836x [3.9758] 




5701 


670 


5.3197 [3.9921] 


lower mantle 


5701 


670 


-6.4761x+5.6283x 2 -3.0807x 3 [4.3807] 




5600 


771 


7.9565 [4.4432] 




5600 


771 


-6.4761x+5.6283x 2 -3.0807x 3 [4.4432] 




3630 


2741 


7.9565 [5.4915] 




3630 


2741 


-6.4761x+5.5283x 2 -3.0807x 3 [5.4915] 




3480 


2891 


7.9565 [5.5665] 


outer core 


3480 


2891 


-1.2638x-3.6426x 2 -5.5281x 3 [9.9035] 




1221.5 


4260.5 


12.5815 [12.1582] 


inner core 


1221.5 


4260.5 


-8.8381x 2 [12.764] 







6371 


13.0885 [13.0885] 



Table 4: Earth density versus the radius; x = t/te, D = ve — r. 



33 




2 4 6 8 10 12 
L (10 3 km) 



Figure 4: Earth density profile. 



For very long baseline experiment with high energy neutrino beams, such as H2B, the matter 
effect will be important. The deepest reach of a neutrino beam is given by 



D = r E {l-^--gr) (73) 

which is 80 km in the case of H2B. The neutrino beam will go through mostly the earth crust but 
also part of the upper mantle. On the average the earth density along the path of H2B varies from 
2.6 g/cm 3 at the initial beam entry into the earth to 3.8 g/cm 3 in the middle point of the beam 
and back to 2.6 g/cm 3 when the beam exits into the detector. Local variation along the beam path 
has to be carefully modeled in order to investigate the CP effect. 



4.1.4 Interaction cross sections 

The detection of the neutrino flavor is through the charge current interaction. For the neutrino 
energy which is small compared to the mass of the W-boson, the charge current cross sections for 
the electron and muon neutrino are given by 

affi = 0.67 x 10~ 38 cm 2 ^(GeV) (74) 

affi = 0.34 x 10- 38 cm 2 ^(GeV) (75) 

For the tau neutrino, the above expression is subject to a threshold suppression. The threshold 
for the production of the tau is 

2 

in 

E T = m T + — ^- = 3.46 GeV (76) 
zrriN 
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The charge current production cross section of the r is usually given numerically as a function 
of the neutrino energy. We fit the numerical cross sections from the threshold to 100 GeV and 
obtained the following expression 



a 



(r) 
uN 



/ a uN 



(E V -E T Y 
c + c\E v + c 2 El 



9(E U — Et), 



(77) 



where cq = —84.988, c\ = 18.317, and C2 = 1.194. As shown in Fig. |g, the fit is good to within 
3%. The difference occurs mostly in the neutrino energy region of 20-40 GeV. The fit is valid for 
E v > 4.0 GeV. 




Figure 5: A fit of the numerical ratio of the r to fi charge current production cross section from the 
threshold to 100 GeV. The boxes are the data and the dotted curve is the fit. The horizontal axis 
is the neutrino energy in GeV and the vertical axis is the ratio of r to fi production cross sections. 



4.1.5 Neutrino Statistics 

The number of neutrino events of E u and flavor /3, from a neutrino beam of energy E u and flavor 
a, to be observed at a baseline L is given by 

N vp = L)u{E v )P a ^{E v , L), (78) 

where <&{E y) L) is the neutrino flux including the detector size, P a ^p{E v , L) is the oscillation prob- 
ability and a(E v ) the neutrino charge current cross section. Without committing to an accelerator 
and a detector design we can not specify the absolute neutrino statistics precisely. In the simulation 
presented in the next subsection, we will take an arbitrary normalization. From the discussion of 
Subsection 4.1.1, we see that the relative energy spectrum of the neutrino beam is specified. For 
the neutrino factory we will not specify the parent muon energy, hence each neutrino energy is 
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taken to be the average energy of the neutrino beam energy as also discussed in Subsection 4.1.1. 
Then we have 



N oc { jjPa^/3(E u ,L) ^ neutrino - factory 

\ (1 - x u ) A (l + be~ Xl,D )^P a ^j3(E v ,L) meson - neutrino 

The dependence of N u . on L and E v depends on the behavior of the oscillation probability on these 
variables. The general behavior of oscillation probability is complicated because it is a function 
of two MSD's which appear in oscillating functions. However, for LBL with neutrino energy of 
order GeV and varying within a relatively small range, the contribution of the solar MSD, Am^, 
is small and proportional to (L/E u ) 2 . The leading contribution comes from the atmospheric MSD 
scale, Ara§ 2 . Hence for L/E v « (Am^) -1 , P a ->p{E v ,L) oc L 2 /E 2 . For L/E u rs (Amf^) -1 , 
P a _>l3(E u , L) is not subject to the (L/E v ) 2 suppression except for the energy dependence entering 
the matter effect which modifies the mixing angles. Therefore, we have 

f E for L/Eu « (Ami)- 1 , , 

\ fir for LjE v » (Am| 2 ) 1 

for neutrino factories and 

f (l-x u ) A (l + 5e-^ D ) for L/E v « (Am^)" 1 
^ I (1 - x u ) A (l + 5e~^ D )§ for « (Am^)- 1 1 j 



for conventional neutrino beams by using Eqs. (|68| , |7l| , |74| , [7q ) . Naively from the above, we expect to 
have better results at higher energy for neutrino factories while at lower energy for conventional 
beams. In the following, we will discuss in detail where are the best place to do measurements for 
some of the most important quantities. 



4.2 Physics of H2B 

In this subsection we will summarize in general terms the physics goals of H2B and compare its 
capability with those of different distances |43]| . A preliminary study of the sensitivity of BAND 
on the measurements of the physics goals, with a specific detector as an example, will be discussed 
in the next section. We assume that in 5 to 8 years the solar neutrino mixing parameters will be 
more accurately determined. Also the parameters of the atmospheric neutrinos will be narrowed 
down. A broad range of physics goals can be defined for H2B, depending on the neutrino beam 
either from a meson source or from a z/-factory. We emphasize the advantages of an experiment 
at a very long baseline, such as 2100 km, and the neutrino energy in the GeV range. To limit the 
scope of the analysis, we will only discuss the scenario of three neutrinos with MSW-LMA as the 



solution to the solar neutrino problem. Some of the following results will appear in Ref. [44]. The 
earth density has been chosen to be a constant of p = 3 g/cm 3 , other mixing parameters are set as 
follows: Am| 2 = Am^ = 3 x lO" 3 ^ 2 , Am^ = 5 x W~ 5 eV 2 , d 12 = 6> 23 = 45°, and 13 = 7°. The 
leptonic CP phase 8 is set to be 90° and the matter effect constant A/E u is 2.3 x 10 eV 2 /GeV. 
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4.2.1 Mixing probability sin 2 26*13 

The oscillation probability P{y^ — ► v e ) is a direct measurement of sin 2 2#i3 since P (y^ — ► z/ e ) oc 
sin 2 26 ) i3. Of course sin 2 2#i3 also appears in f e — > channel, but experimentally it is much more 
difficult. The statistical significance of P{y^ — > f e ) can be measured by the figure of merit defined 
by Piy^ — > u e )/8P(uij, — > v e ), where the error 5P{v^ — * v e ) can be written as 

SP = 5N s /<S>a = 5(N - N b )/$a = ^ 5 2 N + 5 2 N b /$o = ,J N + r 2 N% /<Sht. (82) 

Here N is the total candidate event, N s the signal event and Nf, the background event. The 
uncertainty in estimating the background is represented by a parameter r, which is typically a few 
percent. In this report, we will use r = 0.1. If we express N b = f&a, where f is the background 
fraction in terms of neutrino events of the original flavor, we have 

5N S = ^P§(j + /$<r + r 2 f 2 <S> 2 a 2 . (83) 

It is clear that if P « f, there will be no sensitivity for the given measurement of the signal. 
Generally speaking, at shorter baseline when P is small, the effect of backgrounds is larger. Typical 
values of / for water Cerenkov detectors are a few percent]]. Contributions to / from the beam are 
at the level of 0.6% for meson-neutrino beams. We conservatively choose / = 0.03 for conventional 
beams and / = 0.02 for neutrino factories. 

The figure of merit P/5P can then be written as 

P P P 

6P ~ SNj^ ~ T^T^F 

Using the analytical formula for oscillations in Sec. 3.3.2 and Eqs.(|6S|) and (|74|-|75D, this quantity 
can be plotted as a function of the neutrino energy E at different baselines. 

To see the role played by the background let us first consider the figure of merit by ignoring the 
background, i.e., setting / = 0. The results are given in Fig. ^a and b. As shown in Fig. ^a there 
is no E op t for a given L where the figure of merit is a maximum for E v < 20 GeV for a neutrino 
factory. It would be more advantageous to work at higher neutrino energies and shorter baselines. 
For a conventional neutrino beam, there is an E op t associated to each L, and the smaller baseline 
seemed to have higher figure of merit as shown in Fig. ^3. 

The above picture is changed when the effect of the background is taken into account. As 
shown in Figs, ^c and |6]d, an E op t exists for each distance with the figure of merit increasing with 
the oscillation distance for both the conventional neutrino beams and neutrino factories. For the 
neutrino factory, the figure of merit of the longer baselines such as 2100 km and 3000 km are 
much higher than the shorter baselines of 700 km and 300 km. For the meson-neutrino beam the 
longer baseline is also better although the difference is not so prominent in comparison with the 
case of neutrino factory. This comparison of Figs |6|c and |6|d with Figs ||a and ||b demonstrates the 
importance of the effect of the background in deciding the relative merits at different baselines. 

Since neutrino beams generally have a broad energy distribution, the correct figure of merit 
should integrate over the expected energy spectrum. While at neutrino factories, the neutrino 

1 See next section for the water Cerenkov calorimeter and Ref. S] for water Cerenkov ring imaging detector. 
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Figure 6: Figure of merit for sin 2 26^ measurement at a) neutrino factories without backgrounds, 
b) meson-neutrino beams without backgrounds, c) neutrino factories with f=0.02 and r=0.1, and 
d) meson- neutrino beams with f=0.03 and r=0.1. 
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energy spectrum is well represented by Eq.(p8[), it is more complicated for meson- neutrino beams 
since specific beam design can alter significantly the energy spectrum. As a general rule of thumb, it 
is always better to have a larger width of the figure of merit curve. To this respect, it is advantageous 
at longer baselines as shown in Fig. [| A simple trial effort using Eq.(]68|) for neutrino factories and 
Eq.([7l|) for meson-neutrino beams shows that L=2100 km is preferred. 



4.2.2 Leptonic CP phase 5 

As discussed in the Introduction, the leptonic CP phase for neutrinos may be large as the neutrino 
mixing angles are large. The CP phase 5 can be measured by looking at the difference of the 
oscillation probability between P(i^ — > v e ) and P{pn — * v e ) or between P{y e — > and P{v^ — > 
v e ). While the former has to have the matter effect removed in order to obtain the CP effect, the 
later, under the assumption of CPT conservation, allows the isolation of the matter from the CP 
effect, but it can only be done at a neutrino factory. 
We define 

AP = P(^ -» v e ) - P{vp -^u e ) = P-P (85) 

and the pure CP phase can be obtained by subtracting the matter effect AP(<5) — AP(5 = 0). The 
error of AP can be written as 



6(AP) = vWT^P = JW + <™£ (86) 



and the figure of merit for the leptonic CP phase measurement can be written as 
AP(<5) - AP(5 = 0) _ AP(6) - AP(5 = 0) 



(87) 



using Eq.(||). Fig. |a shows the fi gure of merit for neutrino factories with backgrounds as discussed 
above. It is clear from the figure that, for every L, there is an optimum energy for CP phase 8 
measurement. Once L/E is fixed at the optimum energy, the sensitivity to CP phase 5 is a factor 
of 2 better at L=2100 km than that at L=300 km. The sensitivity to CP for conventional neutrino 
beams is shown in Fig. ^b by using Eq.(|7l"|) and (|74|-[75|). A similar optimum energy also exists, but 
the sensitivity at L=2100 km is about 20% lower than that at L=300 km at their respective peak 
values due to unfavorable beam divergence at long distance. The integrated figure merit, however, 
still favors the longer baseline. 

Similar conclusions can be obtained by using other quantities, such as AP/(P + P), P/P, 
P{y e -> u,j,) - P(z^ -> v e ), etc. 

^,From Eq.(|87|), we can see that the figure of merit is roughly proportional to 1/vP- This means 
that it is much better to measure the CP phase when P is small, i.e., in the v e — Vn or v e — v T 
channel, but not the f^ — y T channel, thanks to the small term sin 2 2#i3 in P(u e — iaJ and P{y e — v T \ 
In reality, we can only use the v e — Vn channel due to higher experimental efficiencies. 



39 




Figure 7: Figure of merit for CP phase measurement at a) neutrino factories with f=0.02 and r=0.1, 
and b) meson-neutrino beams with f=0.03 and r=0.1. 



4.2.3 Sign of Amj 2 

Since the vacuum oscillation probability is an even function of AmL the presence of the matter 
effect is necessary in order to measure the sign of Am| 2 . The simplest way is to compare the 
measured probability Piy^, — ► v e ) with expected P + {v^ — > v e ) and -P~(z^ — ► v e ) where, P + 
assumes Ara^ > and P~ for Am| 2 < 0. The channel 9u — ► P e may be more advantageous to use 
if Am| 2 < 0. This measurement can only be done if sin 2 2#i3 is sizable and the v e or v e appearance 
signal is statistically significant. Other channels, such as the i/„ — * Va survival channel, are less 



sensitive due to the very small matter dependent terms. It has been suggested [45] that total muons 
from all channels at neutrino factories, e.g. the u e (u e ) — ► v^ip^) appearance and v^ipp) — > v^v^) 
survival channels, be used to reduce systematic errors. Here we confine ourselves only to relatively 
simple approaches to illustrate the statistical importance at different baselines. The figure of merit 
in the present case can be written as 

P+ _ p- p+ _ p- p+ _ p- 



6P 5N s /$o x lP±L + r 2 P 



Fig. ^ shows this figure of merit as a function of energy at different baselines for a) neutrino 
factory and b) conventional neutrino beams. It is clear from the plots that in both cases, it is much 
better to have an experiment at L=2100 than at L=300 km. 
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Figure 8: Figure of merit for the sign of Am^ at a) neutrino factories with f=0.02 and r=0.1, and 
b) meson-neutrino beams with f=0.03 and r=0.1. 



4.2.4 Matter effect 

The matter effect can be measured by looking at the difference between Vn — > v e and Vn — » P e , 
although non-zero CP phase 5 can alter the result by as much as 20%. The survival channels 
Vy. —* v a and — > are not very effective for the present consideration since matter effects in 
these channels appear only in the non-leading terms. It is better to use the channels that involve 
v e or v e . 

Denote P\ = P{v^ — > v e ) and P 2 = P(vu — * P e )> the relevant figure of merit can be written as 



P1-P2 
5{P l - P 2 ) 



P1-P2 



Px-Pi 



P1-P2 



V<5 2 Pi + 5 2 P 2 



S 2 N U 



5 2 Np 



Pi+f 1 r 2 f 2 1 p 2+f 1 r 2 f 2 



Fig. ^ shows this figure of merit as a function of energy at different baselines for a) neutrino factory 
and b) conventional neutrino beams. It is clear from the plots that in both cases, it is much better 
to have an experiment at L=2100 than at L=300 km. 



4.2.5 Precision measurement of Am| and sin 2 29 



7 28 



Although A77132 and sin 2 2#23 have been measured at super-K and hopefully will be confirmed by 
K2K, MINOS, OPERA and ICARUS, it is still interesting to measure them at different baselines 
and possibly improve the precision. Am 2 2 and sin 2 2#23 can be directly related to the survival 
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Figure 9: Figure of merit for the matter effect at a) neutrino factories with f=0.02 and r=0.1, and 
b) meson- neutrino beams with f=0.03 and r=0.1. 



probability Piy^ — * v^). The figure of merit can be written as 

1 -P 1 - P 1-P 



SP 5N s /*a ./m + r2/2 



(90) 



Fig. 1C shows the figure of merit at different baseline for a) neutrino factories and b) conventional 
neutrino beams. Since the background for muon identification is generally much smaller than in 
the case of the electron, we use f=0.01 for both neutrino factories and conventional beams. A nice 
peak corresponding to Am^I where the oscillation is at a maximum can be found and its position 
is a good measure of Amjj 2 . It is clear that for neutrino factories it is better at 2100 km and for 
conventional beams, it is better at 300 km. The integrated figure of merit at L = 2100 km is much 
higher than that at L = 300 km for a neutrino factory. For the meson-neutrino beam the longer 
baseline is also higher. 

4.2.6 f M — > v T Appearance 

The appearance of v r from a Vn beam is an unambiguous proof of the — > v T oscillation. Precision 
measurement of this oscillation probability is crucial in establishing the oscillation patterns and in 
determining whether or not Va — ► v T is the only dominant oscillation mode or there are still rooms 
for the — > v s oscillation. Although we expect that in the next few years, K2K and/or OPERA 
will observe the production of r because of the large Vn — > v T probability, it is desirable that BAND 
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Figure 10: Figure of merit for Am| and sin 2 2#23 at a) neutrino factories with f=0.01 and r=0.1, 
and b) meson-neutrino beams with f=0.01 and r=0.1. 



possesses a good r identification capability. We assume that our detector can identify the r based 
on statistics methods, namely event selection via signature on kinematics. The figure of merit can 
be written as 



V T ) 



p 



6N s /$a 



P+f 

<E>cr 



+ r 2 / 2 



(91) 



Fig. 11 shows this figure of merit as a function of energy at different baselines for a) neutrino 
factories and b) conventional neutrino beams. Although the neutrino beam does not contribute to 
the background, it is expected to be high after event selection. Hence we set f=0.03 for both neutrino 
factories and meson-neutrino beams. Due to the threshold of tau production, the experiment has 
to be done at E„ > 4GeV. That is different from what is needed for CP phase measurements at 
L < 2000 km. While shorter baseline seems better at a given energy, it is generally better to run 
only one energy for all measurements to accumulate enough statistics. To this end, L=2100 km is 
much better than that at L=300 km. 



4.2.7 Exotic scenarios 

• Sterile neutrinos Although the Super-K collaboration has demonstrated that the sterile 
neutrino is excluded at the 99% CL to be the dominant oscillation partner of v^, it has not 
ruled out a sizable oscillation probability for — ► v s |4(| . The presence of — ► v s oscillation 
will lead to a suppressed survival probability v „ — > i/„ at E ~ 15 GeV due to matter effects. 
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Figure 11: Figure of merit for tau appearance at a) neutrino factories with f=0.03 and r=0.1, and 
b) meson-neutrino beams with f=0.03 and r=0.1. 



Then very long baseline experiments, such as H2B, will have the advantages to search for the 
suppression. 

• Neutrino decay The possibility of neutrino decay as an explanation of Super-K's atmo- 
spheric neutrino observations has not been completely excluded p6| . However, it has been 
predicted that there is a significant difference on the survival probability P{ya — > Vn) 
between the oscillation and decay scenarios at L/E u (km./GeV) ~ 400. Therefore the figure 
of merit to make the discrimination of the two scenarios is equivalent to that given in Fig. |io| . 
Hence we expect that K2K and MINOS will be able to answer the neutrino decay question 
in the near future before we are online. 



4.2.8 Summary 

In summary, we found that by including the effect of backgrounds, the figure of merit for the very 
long baseline of L = 2100 km is higher than that at L = 300 km for all the measurements that we 
have considered for a neutrino factory. For the meson-neutrino beam L = 2100 km is still higher 
in most of the measurement. Hence according to the figures of merit, it is clearly advantageous to 
do experiments at a longer baseline such as L = 2100 km. In Table || we list all the figures of merit 
obtained from the maxima of the various curves given in Figs. [7-11, 



To put things in perspective, a few remarks are due. A 3a statistical significance corresponds 
to value 3 in this table but unfortunately, the values of the figures of merit as given in Table [| 
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neutrino factory 


meson-neutrino beam 




300 km 


2100 km 


300 km 


2100 km 


sin^ 26*13 


3.5 


10.5 


6.0 


8.0 


CP phase S 


0.5 


0.95 


0.62 


0.5 


sign of Am 2 2 


1.0 


10. 


2.5 


8.0 


matter effects 


0.2 


8. 


0.5 


6.0 


Am^ and sin 2 2# 2 3 


180 


450 


500 


350 


r appearance 


10 


50 


9.5 


13.5 



Table 5: Summary of relative figure of merits for various measurement at baselines L=300 and 2100 
km. 



are only relative. Eventually the absolute statistics should be used to obtain the absolute figure 
of merit. There are two absolute normalization factors which we treat as free parameters in the 
present discussion: one for neutrino factory and the other for conventional neutrino beams. These 
parameters depend on the beam intensity, the detector mass, efficiency, etc... In addition, we also 
need to know accurately the oscillation parameters such as sin 2 (2# 13 ), S, Am 2 2 , etc. However 
independent of these parameters it is clear from the table that at neutrino factories, L=2100 km 
is always better than L=300 km. For the meson-neutrino beams, it is better at L=2100 km for 
sin 2 26i3, the sign of Am|, the matter effects and tau appearance. For CP phase measurement, 
it is comparable at all different baselines while for Am^ and sin 2 (26*23), L=300 km has a higher 
figure of merit. However, when the integrated figures of merit are considered the L = 2100 is always 
preferred. 
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5 Far detector 



As we discussed before, our main physics objectives include the measurements of sin #13, Am]], 
the leptonic CP phase 5 and the sign of Am^ . All of these quantities can be obtained through the 
survival probability P(f M — > v^) and the appearance probability P(f A1 (f e ) — > f e (f M )) and P(P A1 (i/ e ) — > 
Veipy)). To measure these quantities, a detector should: 1) be able to identify charge leptons: e, 
and r; 2) have good pattern recognition capabilities for background rejection; 3) have good energy 
resolution for event selection and for determining P Q ,^ ( g(E) as a function of the neutrino energy; 4) 
be able to measure the charge for in the case of v factories; and 5) has a large mass (100-1000 kt) 
at an affordable price. 





Iron 


Liquid 


Water Ring 


Under Water /Ice 




Calorimeter 


Ar TPC 


Imaging 


Cerenkov counter 


Mass 


5-50 kt 


1-10 kt 


50-1000 kt 


100 Mt 


Charge ID 


Yes 


Yes 


? 


No 


E resolution 


good 


very good 


very good 


poor 


Examples 


Minos 


ICARUS 


Super-K, Uno 


Amanda, Icecube 




Monolith 




Aqua-rich 


Nestor,Antares 


Price ($) 


lOM/kt 


30M/kt 


2M/kt 


lK/kt 



Table 6: Currently proposed detector for v factories and conventional v beams. 



Currently there are four types of detectors proposed [47, ^1, as listed in table ||, for large LBL 
experiments. As seen from the table, they are generally expensive and some of them are difficult 
to have magnet in order to measure the sign of the electric charge of the muon. However, much 
cheaper detectors, say at a price target of 0.2 M/kt, are not impossible. They are very attractive 
since they allow us to build a detector of 100 - 1000 kt without much difficulties. Here we describe 
one such option which was first proposed in ref. [49]. Of course a significant amount of R& D are 
needed for this option and we are open to other detector designs. 



5.1 Water Cerenkov Calorimeter 

Water Cerenkov ring image detectors have been successfully employed in large scale, for obvious 
economic reasons, by the 1MB and the Super-Kamiokande collaborations. However a substantial 
growth in size beyond these detectors appears problematic because of the cost of excavation and 
photon detection. To overcome these problems, we propose here a water Cerenkov calorimeter with 
a modular structure, as shown in Fig. |l2[ 

Each tank has dimensions 1 x 1 x 10m 3 , holding a total of 10 t of water. The exact segmen- 
tation of water tanks is to be optimized based on the neutrino beam energy, the experimental 
hall configuration, the cost, etc. For simplicity, we discuss in the following a module of 1 m in 
thickness, corresponding to 2.77 Xo and 1.5 Ao- The water tank is made of PVC with aluminum 
lining. Cerenkov light is reflected by the aluminum lining and transported towards the two ends of 
the tank, which are covered by wavelength shifter(WLS) plates. Light from the WLS is guided to 
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Figure 12: Schematic of water Cerenkov calorimeter 




Figure 13: Schematic of a water tank. 
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a 5" photon- multiplier tube(PMT), as shown in Fig. 13. The modular structure of such a detector 
allows it to be placed at a shallow depth in a cavern of any shape(or possibly even at surface), 
therefore reducing the excavation cost, in particular for the presently planned H2B site where a hall 
is available and a large scale excavation is unnecessary. The photon collection area is also reduced 
dramatically, making it possible to build a large detector at a moderate cost, which is targeted at 
$0.2M/kt. 

A through-going charged particle emits about 20,000 Cerenkov photons per meter. Assuming a 
light attenuation length in water of 20m and a reflection coefficient of the Aluminum lining of 90%, 
we obtain a light collection efficiency of about 20%. Combined with the quantum efficiency of the 
PMT(20%), the WLS collection efficiency(25%) and an additional safety factor of 50%, the total 
light collection efficiency is about 0.5%. This corresponds to 100 photoelectrons per meter, which 
can be translated to a resolution of 4.5%/\/E. This is slightly worse than the Super-Kamiokande 



detector and liquid Argon TPC but much better than iron calorimeters [47]. 

If this detector is built for a v factory, a tracking device, such as Resistive Plate Chambers 
(RPC)j50| will be needed between water tanks to identify the sign of charge. RPCs can also be 
helpful for pattern recognition, to determine precisely muon directions, and to identify cosmic- 
muons for either veto or calibration. The RPC strips will run in both X- and Y-directions with a 
width of 4 cm. A total of ~ 10 5 m 2 is needed for a 100 kt detector, which is more than an order 
of magnitude larger than the current scale[50|. R&D efforts would be needed to reduce costs. 



The magnet system for such a detector can be segmented in order to minimize dead materials 
between water tanks. If the desired minimum muon momentum is 5 GeV/c, the magnet must be 
segmented every 20 m. Detailed magnet design still needs to be worked out; here we just present a 
preliminary idea to initiate the discussion. A toroid magnet similar to that of Minos, as shown in 
Fig. m can produce a magnetic field B > 1.5 T, for a current I > 10 4 A. The thickness of the magnet 
needed is determined by the error from the multiple scattering: AP/P = 0.0136\/X/Xo/0.3BL, 
where L is the thickness of magnet. For L=50 cm, we obtain an error of 32%. The measurement 
error is given by AP/P ~ 5a/ a = <rP/0.3rBL, where r is the track length before or after the magnet 
and a is the pitch size of the RPC. For P=5 GeV/c, a = 4 cm and r=10 m, the measurement error 
is 9%, much smaller than that from multiple scattering. It should be noted that P^ is also measured 
from the range. By requiring that both P^ measurements are consistent, we can eliminate most of 
the fake wrong sign muons. The iron needed for such a magnet is about 20% of the total mass of 
the water. 

The cost of such a detector is moderate compared to other types of detectors, enabling us to 
build a detector as large as 100 - 1000 kt. The combination of size, excellent energy resolution and 
pattern recognition capabilities makes this detector very attractive. 

5.2 Performance of Water Cerenkov Calorimeter 

A full GEANT Monte Carlo simulation program and the Minos neutrino event generator have been 
used to study the performance of such a detector. Fig. |l^ shows a typical charge current (CC) 
event in our detector. 



A CC v signal event is identified by its accompanying lepton, reconstructed as a jet. Fig. [16 
shows the jet energy normalized by the energy of the lepton. It can be seen from the plot that 
leptons from CC events can indeed be identified and the jet reconstruction algorithm works properly. 
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RPC 

Figure 14: Schematic of a toroid magnet 



Figure 15: A typical Up CC events in the detector. 



49 



It is also shown in the figure that the energy resolution of the neutrino CC events is about 10% in 
both cases. 




to/E v 



Figure 16: The reconstructed jet energy and total visible energy. The fact that Ej et /Ei epton peaks 
around one shows that the jet reconstruction algorithm finds the lepton from CC events. The 
fraction of total visible energy to the neutrino energy indicates that we have an energy resolution 
better than 10% for all neutrinos. The bias is due to invisible neutral hadrons and charged particles 
below Cerenkov thresholds. 



The meson-neutrino beam profile from HIPA has been discussed in Sec. 4 and the energy 
spectra of visible v„ CC events are shown in Fig. 17, A total of at least 30 event/kt per year is 
expected at BAND. For the present report we will not discuss the detector performance at the 
neutrino factory. 

The neutrino CC events are identified by the following 5 variables: E max /Ej e t, L s hower/Ej c t, 
Ntank/Ejet, Rx y /E to t, and R™ ax /E tot , where E jet is the jet energy, E tot the total visible energy, 
E m ax the maximum energy in a cell, L s h OW er the longitudinal length of the jet, N tan k the number of 
cells with energy more than 10 MeV, R xy the transverse event size and R^* the transverse event 



size at the shower maxima. Fig. |18| - |22| shows these variables for all different neutrino flavors. It can 



be seen that v e CC events can be selected with reasonable efficiency and moderate backgrounds. 
Table shows the final results from this pilot Monte Carlo study. For v e and events, v T 
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Figure 17: Beam profile of HIPA-Beijing 
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Figure 18: The maximum energy deposit in a cell normalized to the jet energy, E max /Ej ct , for 
various type of neutrino events. 
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Figure 19: The longitudinal length of a jet normalized to the jet energy, L s h wcr/Ej e t, for various 
type of neutrino events. 
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Figure 20: The number of cells fired normalized to the jet energy, N ta nk/Ej c t, for various type of 
neutrino events 
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Figure 21: The transverse event size normalized to the total visible energy in the detector, R xy /E tot , 
for various type of neutrino events. 
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Figure 22: The transverse event size at the shower maxima normalized to the total visible energy 
in the detector, R^ ax /E to t for various type of neutrino events. The distribution of u e is different 
from all the others. 
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CC events are dominant backgrounds, while for v T , the main background is v e . It is interesting 
to see that this detector can identify v T in a statistical way. Similar results are obtained for a 
detector with 0.5m water tanks without RPCs. These results are similar to or better than those 
from water Cerenkov image detectors [28] and iron calorimeters [51]. we would like to point out that 
using sophisticated jet reconstruction algorithm, shower shape analysis and neural network, better 
results are expected. 









v T 


CC Eff. 


30% 


53% 


9.3% 


v e CC 




>1300:1 


3:1 


v e NC 


166:1 


665:1 


60:1 


^ cc 


700:1 




270:1 


v„ NC 


92:1 


> 6000:1 


39:1 


u T CC 


20:1 


12:1 




v T NC 


205:1 


1100:1 


61:1 



Table 7: Results from Monte Carlo simulation: Efficiency vs background rejection power for differ- 
ent flavors. 



Using Table [7|, we can explore the sensitivity of our detector to various quantities, such as 
sin 2 2#i3, CP phase term, etc. Fig. 23 shows the sensitivity to sin 2 2#i3 for 2 years running of a 
100 kt detector. If systematic errors can be controlled so that it does not dominate, we can reach 
sin 2 2<9i3 ~ 0.006, corresponding to 13 ~ 2.2°. 



5.3 Other physics opportunities 

The detector we proposed above has an extreamly large mass, very good energy and angular 
resolution, hence a rich physics programs other than long baseline neutrino oscillation experiments 
can be studied. We discuss briefly other physics opportunities in the following, obviously more 
detailed studies are needed. Depending on the over burden of the detector, some or all of the 
following incomplete list of physics topics can be explored: 

• atmospheric neutrinos A large data sample of atmospheric neutrinos with very good energy 
and angular resolution can be collected. We expect to improve the results obtained by Super- 
K. 

• Neutral current cross section The charge current interaction cross section is well-know, 
but there is very little data on neutral current cross sections. The BAND detector can 
contribute this area. 

• supernova If we can push the detection threshold down to about 10-15 MeV, this detector 
can observe supernovae at distances from us up to hundreds of kpc. 

• primary cosmic-ray composition By measuring multiple muons, the primary cosmic-ray 
composition, dominated either by heavy nuclei or light nuclei, can be determined. 
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Figure 23: Sensitivity to sin 2 26 13 at H2B. 
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UHE cosmic rays With an over burden of a few hundreds meters of water equivalent, 
muons with energy above 50 GeV can be detected. This is an ideal region for studying the 
composition of primary cosmic rays beyond the "knee" and the details of hadronic interactions 
in the fragmentation region. 

searches for dark matter WIMP's as cold dark matter can annihilate in the core of the 
sun and the earth. By looking at excess of muons from the core of the sun and the earth, we 
can search for signs of dark matter. For a few 100 kt detector, the allowed region from the 
DAMA experiment can be accessed. 

searches for monopoles Monopoles are predicted by many theories beyond the Standard 
Model. It can be searched for by looking at slow moving particles with high dE/dx in our 
detector. A few orders of magnitude improvement over results from MACRO are expected, 
searches for point sources High energy neutrinos and muons from point sources other 
than galaxies can provide important information about the evolution of the universe and 
the acceleration mechanism of cosmic-rays. Our detector can contribute significantly in this 
respect. 

searches for exotic particles The history of physics is full of consummated surprises. 
"Unexpected" and "exotic" particles may be detected in our detector, such as fractionally 
charged particles. 
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